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We calculate perturbative Wilson loops of various sizes up to loop order n = 20 at different lattice 
sizes for pure plaquette and tree-level improved Symanzik gauge theories using the technique of 
Numerical Stochastic Perturbation Theory. This allows us to investigate the behavior of the per- 
turbative series at high orders. We observe differences in the behavior of perturbative coefficients 
as a function of the loop order. Up to n = 20 we do not see evidence for the often assumed facto- 
rial growth of the coefficients. Based on the observed behavior we sum this series in a model with 
hypergeometric functions. Alternatively we estimate the series in boosted perturbation theory. 
Subtracting the estimated perturbative series for the average plaquette from the non-perturbative 
Monte Carlo result we estimate the gluon condensate. 
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1. Introduction 

Since the introduction of the non-perturbative gluon condensate by Shifman, Vainshtein and 
Zakharov [jl|] there have been many attempts to obtain reliable numerical results for this quantity. 
Soon it became clear that lattice gauge theory provides a promising tool to calculate it from Wilson 
loops Wnm of sizes N xM.ln [§, ^|] the plaquette was used whereas larger Wilson loops have been 
investigated in [|], |5||. In all cases it turned out to be crucial to know the perturbative tail of the 
Wilson loops as precisely as possible. In the last decade the application of Numerical Stochastic 
Perturbation Theory (NSPT) ^ pushed the perturbative order of the plaquette up to order n = 10 
[§] and even n = 16 [@|. 

Apart from the extraction of the gluon condensate there is a general interest in the behavior 
of perturbative series in QCD (for a recent investigation see It is generally believed that 
these series are asymptotic, and assumed that for large n the leading growth of the coefficients a„ 
is factorial. Using the technique of NSPT one reaches orders of the perturbative series where a 
possible set-in of this assumed behavior can be tested - at least for finite lattices. In a more recent 
paper JTo| ] Narison and Zakharov discussed the difference between short and long perturbative 
series and its impact on the determination of the gluon condensate. 

In order to test to what extent &{a 2 ) improvement influences the behavior of the series we 
used in addition to the standard Wilson plaquette gauge action the tree-level improved Symanzik 
gauge action obtained by Weisz implementing Symanzik's improvement programme [juj]. 

We present perturbative calculations in NSPT up to order n = 20 for Wilson loops using the 
Wilson action for lattice sizes L 4 with L = 4,6,8,12. In case of the Symanzik action we have 
computed the W^m f° r L = 4, 6, 8, 10. The computation for L = 12 were performed on a NEC SX-9 
computer of RCNP at Osaka University, for L = 10 at the HLRN Berlin/Hannover, all others on 
Linux/HP - clusters at Leipzig University. 



2. NSPT calculation up to n = 20 

NSPT allows perturbative calculations on a finite lattice up to loop order n which practically 
cannot be reached by the standard diagrammatic approach. A limit is set only by storage limitations 
and machine precision. The algorithm is introduced and discussed in detail in [Q, 12] - we will not 



present it in this paper. For a detailed discussion of our results we refer to a forthcoming paper [jOy. 

In order to fix the notation we write the general expansion of a Wilson loop of size N xM in 
terms of the bare lattice coupling g as 

«=0 



The coefficients W NM are determined with NSPT and need to be known to a very good precision. 
In Figure [l] we show some results for these coefficients using the Wilson and Symanzik actions. 
One can see from that Figure that the statistical errors within the NSPT algorithm are rather small. 

At all n all the loops behave rather similarly, with the coefficents decreasing smoothly, in a 
similar way for all loop sizes. At smaller n we see that the coefficients change sign and magni- 
tude in a rather unpredictable way, particularly for large loops. We recognize some remarkable 
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Figure 1: Selected for L — 12 and Wilson action (left) and for L = 10 and Symanzik action (right). 
Positive/negative signs of the coefficients are given by open/full symbols. 



differences in the behavior between the two actions. First, for the Symanzik action the sign change 
in the expansion coefficients sets in at smaller Wilson loop sizes N x M and is extended to larger 
loop orders n. Second, the values of the Symanzik coefficients themselves are much smaller than 
their Wilson action counterparts. However, concerning the relative convergence behavior of the 



perturbative series (2A) one should have in mind that the coupling g 2 used in computing physi- 
cal observables is larger in the Symanzik case than in the Wilson case. At larger n a seemingly 
asymptotic behavior (without sign changes) emerges. 

One essential test for the validity of our perturbative calculations with NSPT is the signal/noise 
ratio. The expansion of the (integer) loop-order n is constructed from even powers g 2n . Likewise 
one can build non-loop contributions which are half-integer and combined from odd powers of g - 
they should vanish. In Fig. ^| we show that the loop contributions are always clearly separated in 
magnitude from the non-loop (noise) terms. 

In addition we have to raise the question about the infinite volume limit L — > 00 of the loop 
series. This is related to the definition of the gluon condensate which is clearly a L = 00 quantity. 
We follow the approach to extrapolate each coefficient at finite L from wjf^ L to wj^ „ using the 
ansatz 



W, 



(») 



NM.L 



®NM,n + bf/M,n L 4 + CfjM,n log(L) L 



(2.2) 



3. Perturbative series at large order 



It is generally assumed that perturbative series in continuum QCD are asymptotic. The sit- 
uation might be different for perturbative series on finite lattices. Here we have both ultraviolet 
and infrared cut-offs and the series could be summed up to a finite value. With our computed 
coefficients up to order n = 20 we are able to check this conjecture to a so far unrivaled level. 

To perform such a summation we propose a model which is a generalization of the ansatz used 



in [14]. Denoting the generic coefficients of the perturbative series as c n we have found that a large 
set of our data above some loop order «o can be described rather well by the following ansatz for 
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Figure 2: Selected W^" versus n together with typical sizes (in magnitude) of non-loop contributions. 
Positive/negative signs of the coefficients are given by open/full symbols. Left: L = 12 and Wilson gauge 
action, right: L = 10 and Symanzik action. Points at half-integer n (the non-loop terms) are pure noise - they 
show a good signal/noise ratio. 



the ratio of subsequent coefficients 



«.l-l±?l+ 



n(n + s) ' 



n > no , 



(3-1) 



where u,y,p,s are free parameters to be fitted. The series based on (3.1) can be summed up to 
n = oo for g 2 < 1 / u. For no = 1 we obtain as solution 



(s+l)W^ 2 F l ((l+s-Y)/2-T,(l+s-Y)/2 + r,s+l;g 2 u) 

p-(s + l)uj 



Wnm,< 



T = ^{{Y+s+l) 2 u-Ap)/u. 
As final hypergeometric model at no = 1 we use 



(3.2) 



(1) 

W NM ,HYP = W NM ,oo - (citfYP ~ WfjM)g 



(3.3) 



where c\^hyp is the coefficient of g 2 in the corresponding g 2 -expansion of Wnm,°°- For other values 



of «o the formulae (3.2) and (3.3) change accordingly. In Fig. H we show two example plots for 



the ratio r„ together with the model fit based on (3.1 ). In all investigated cases we do not observe a 



factorial growth of the coefficients up to order n = 20. Formally, one can use model ( J3-3| ) (used for 
finite lattices) also for the L — > oo extrapolated coefficients. 

A possible alternative method consists in applying boosted perturbation theory, i.e. a rear- 
rangement of the series in terms of a boosted coupling gt, [||] (Wu = P)- Using the summed 
perturbative series P n *(g 2 ) = 1 + L"=i s" we § et the boosted summed plaquette Pio.bigb) as 
follows 



■8b 



20 



P 2 0^) = l + K I ) 1 gi" 



P20(g) 



(3.4) 



B=l 
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Figure 3: Domb-Sykes plots for W\\. Left: L = 12 and Wilson action, right: L = 8 and Symanzik action. 



In Fig. |] we compare the summed perturbative series P n *(g) with P n *,b{gb) as function of maximal 
loop order n* for L = 12. One clearly recognizes a distinct plateau reached by boosting for n* < 20. 



The naive series is far from the model value computed from (3.3). Although this behavior of the 
boosted perturbation theory is found for finite L, we assume that it remains valid also in the infinite 
volume limit. 
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Figure 4: Perturbative plaquette P n * as function of loop order n* for naive and boosted series; left: all orders, 
right: zoomed in by the choice 8 < n* < 20. The choice g 2 = 1 .04 corresponds to the largest coupling where 
the solution (3.2) remains real. 

4. Non-perturbative gluon condensate 

Having under control the large order perturbative series we estimate now the non-perturbative 
gluon condensate ((a s /n)GG). Being a physical quantity we should use the L = oo extrapolated 
values of the perturbative coefficients. On the other hand it could be useful to get numbers for the 
"gluon condensate" at finite L also: they provide some numerical insight how the infinite volume 
limit is approached. 

We use the relation between the plaquette measured in Monte Carlo simulations Pmc and its 
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perturbative analogue P: 



(4.1) 



which defines the gluon condensate. In ( |Q| ) /3(g) denotes the standard /3 -function with & being 
its leading coefficient. From ( fi~l| ) it is obvious that the precision of the difference AP = P — Pmc 
determines the reliability of an estimate for (^GG). 



It has been proposed [10] to extract the gluon condensate taking into account a perturbative 



series of the plaquette summed to an order n* (P n *) from the following ansatz 



if 



-bog 2 



(-GG) 

71 



36 V P(g) 

where AP n * contains possible a 2 and a 4 contributions 



AP n *=P„*-Pi 



MC; 



AP„ 



C2{n*) a 2 + cJn*) a 4 . 



(4.2) 



(4.3) 



The authors in Q10| ] argued that non-zero values of the coefficient C2(n*) w& an artefact due to a 
truncation of the series: above some value of n* this coefficient should vanish. Our perturbative 
series for both, the naive series (finite L) and boosted series (L = oo) confirm their hypothesis: for 
boosted perturbation theory the a^-term decreases rapidly to zero for n* > 10 (see left of Fig. ||); 
for naive perturbation theory and finite L one has to sum the series to much larger n*. 
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Figure 5: Left: Comparison of C4(n*) and C2(n*) for naive and boosted (b) perturbation theory in the 
extrapolated limit L = oo. Right: AP/ (a/ro) 2 as function of (a/ro) 2 for boosted perturbation theory. The 
line is a linear fit in a 2 for the range 0.02 < (a/ro) 2 < 0.08, ro denotes the Sommer scale. 

From the results obtained so far we conclude that the most reliable estimate for the gluon 
condensate can be obtained for the Wilson action using boosted perturbation theory. In the right 
of Figure || we show AP/(a/ro) 2 as function of (a/ro) 2 using boosted perturbation theory with 
n* = 20 and our own non-perturbative plaquette measurements at different /3 values and lattice 
sizes for ro/a > 2.5. We clearly recognize a dominant linear a 2 behavior. The necessary conversion 
/3 <->• (a/ro) is accomplished for the plaquette action by the Necco-Sommer relation JIB}]. 

In Table [j] we give numbers for the quantity ((oc. v / 7l)GG) for finite L and boosted perturbation 
theory at L = oo. For the perturbative values at finite L we use the hypergeometric model as given 
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in (3.3). The table shows that the boosted result for extrapolation L = oo agrees within errors with 
the model value at finite lattice size L = 12. The given errors have been estimated by varying the 
fit window j3 m j n < jS < /3 max . Our reference value for the gluon condensate is that for L = oo and 



L 


4{%GG) 


(f GG) [GeV 4 ] 


fit range 


6 
8 

12 


1.08(10) 
1.22(11) 
1.29(17) 


0.034(3) 
0.039(4) 
0.041(6) 


5.73 < j8 < 6.00 
5.78 </3< 6.10 
5.78 < j3 < 6.27 


oo, boosted 


1.33(7) 


0.042(2) 


5.78 </3 < 6.17 



Table 1: Gluon condensate at finite L for Wilson gauge action (ro = 0.467 fm). The last column shows the 
used j3-nt range. 

boosted perturbation theory given in Table [j]. With the chosen ro it is larger than the value obtained 
by sum rule calculations M\. 
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